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ABSTRACT 

We report on an X-ray observation of the Be X-ray Binary Pulsar RX J0059.2-7138, per¬ 
formed by XMM-Newton in March 2014. The 19 ks long observation was carried out 
about three months after the discovery of the latest outburst from this Small Magellanic 
Cloud transient, when the source luminosity was Lx erg s“^. A spin period of 

Pspin=2.762383(5) s was derived, corresponding to an average spin-up of Apin = — (1-27 ± 
0.01) X 10“^^ s from the only previous period measurement, obtained more than 20 years 
earlier. The time-averaged continuum spectrum (0.2-12 keV) consisted of a hard power-law 
(photon index ^^0.44) with an exponential cut-off at a phase-dependent energy (^20-50 keV) 
plus a significant soft excess below ^0.5 keV. In addition, several features were observed in 
the spectrum: an emission line at 6.6 keV from highly ionized iron, a broad feature at 0.9-1 
keV likely due to a blend of Fe L-shell lines, and narrow emission and absorption lines consis¬ 
tent with transitions in highly ionized oxygen, nitrogen and iron visible in the high resolution 
RGS data (0.4-2.1 keV). Given the different ionization stages of the narrow line components, 
indicative of photoionization from the luminous X-ray pulsar, we argue that the soft excess 
in RX J0059.2-7138 is produced by reprocessing of the pulsar emission in the inner regions 
of the accretion disc. 

Key words: accretion - stars: neutron - X-rays: binaries - X-rays: individual (RX J0059.2- 
7138) 


1 INTRODUCTION 

The dominant component in the spectra of X-ray binary pulsars 
is usually well described by a hard power-law model with an ex¬ 
ponential cut-off, sometimes with Fe Ko emission lines. However, 
at low energies several sources show a significant excess over the 
mai n power-law, which has been desc ribed with different models 
(see lLa Palombara & Mereghetti 200^ . for an overview). Although 
various physical pr ocesses have been in voked to produce it, its ori¬ 
gin is still unclear. iHickox et alj ^2004) showed that the origin of 
the soft excess is related to the source total luminosity and that a 
spectral component of this type has been detected in all the X-ray 
pulsars with a sufficiently high flux and small absorption. In fact, 
most of the soft excess sources are at small distances and/or away 
from the Galactic plane. This suggests that the presence of a soft 
spectral component could be a very common, if not an ubiquitous, 
feature intrinsic to X-ray pulsars. Since the study of the soft part of 
the spectrum in Galactic sources is affected by the interstellar ab- 
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sorption present in the Galactic plane, the best sources to study this 
feature are the transient pulsars in the Magellanic Clouds (MCs). 
They can reach high luminosities {Lx > 10^® erg s“^), are at well 
known distances, and, thanks to the low absorption in the MCs di¬ 
rection, can provide high-statistics spectra at low energies. In par¬ 
ticular, the Small Magellanic Cloud (SMC) hosts more than 100 
Be/X-ray Binary systems, about 70 of which show periodic pul- 
sations (see, e.g. iHaberl et alj|201 j : ISturm et'^l2013l : IkIus et aD 
|2014|) . 

One of the most interesting pulsars in the SMC is RX J0059.2- 
7138, also known as SX P 2.76, which was discovered in 1993 by 
ROSAT iHughes|[l994h . A Be type star wa s found in its error re¬ 
gion a nd proposed as a candidate counterpart iSouthwell & Charles! 
Il996h . Optical photometry of this star revealed a periodicity 
of 82 days, likely related to the orb ital period of the system 
jSchmidte et alj2006LlBird et al.ll2012h . The ROSAT X-ray spec¬ 
trum of RX J0059.2-7138 was unusually soft, consisting of a low 
temperature black-body component {kT = 35 eV) plus a steep 
power-law (photon index r=2.4). The presence of a soft excess 
was confirmed by a simultaneous ASCA observation, which also 
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Swift/BAT (15-50 keV) 



Figure 1. RX J0059.2-7138 outburst as observed by Swiftf^K^ (15-50 
keV). The bin time is one day. The BAT light curve is publ icly available 
at ht tp://SWIFT.GSFC.NASA.GOV/RESULTS/TRANSIENTS/ jKrimm et al] 
[5^. The date of our XMM-Newton observation is indicated by the 
vertical dashed line. 


detec ted X-ray pulsations with a period of 2.76 s dKohno et al.l 
l200(]h . The simultaneous analysis of the ROSAT and ASCA data 
showed that the soft component could be described with either 
emission from a thin thermal plasma (kT=0.37 keV) or a broken 
power-law, while a hard power-law (r=0.43) described well the 
spectrum above 2 keV. Both components had a phase-averaged lu¬ 
minosity of erg s“^ (0.1-10 keV): the hard component dom¬ 
inated above 2 keV and showed sinusoidal pulsations with a pulsed 
fraction of ~37%, while no pulsations we re detected in the soft 
component, which dominated below 1 keV dKohno et al.ll200(ih . 

After the 1993 observations, no other bright outbursts were 
seen from RXJ0059.2-7138 until 2013, except for two detections 
of the source with P C A RX TE in Septem ber 2002 and in Novem¬ 
ber 2004 reported bv lOalache et all d2:008l) . The 2013 outburst was 
discovered with the Swift/BAT instrument , which detected th e 
source starting from December 26 (see Fig. m lKrimm et alJl20T4ll . 
A couple of weeks later, INTEGRAL monitorin g observations o f 
the SMC confirmed the bright status of the source i Coe et al.l2014h . 
After the discovery of this outburst, we triggered our XMM- 
Newton Target of Opportunity program aimed at investigating 
the soft excess in high luminosity Be/X-ray transients (hereafter 
Be/XRTs) thus obtaining an observation of RX J0059.2-7138 in 
March 2014. Here we report the results of this observation. 


2 OBSERVATIONS AND DATA REDUCTION 


RX J0059.2-7138 was observed by XMM-Newton between 
2014 March 14 and 15, for a net exposure of ~19 ks (see Table[T] 
for the observation log). In Fig. [T] we show the time of the XMM- 
Newton observation in the context of the outburst as observed by 


Swift/BAT (15-50 keV). 

The XMM-Newton Observatory carries three 1500 cm^ 
X-ray telescopes, each with an European Photon Imaging Cam¬ 
era (EPIC; 0.2-12 ke V) at the fo cus. Two of the EPIC detec¬ 
tors use MOS CCDs dTumer et al.ll200l1) and one uses pn CCDs 
dStruderet al.ll2001). Reflection Grating Spectrometer (RGS; 0.4- 
2.1 keV) arrays dden Herder et all [200 ih are located in the tele¬ 
scopes with MOS cameras at their primary focus. The EPIC cam¬ 


eras were operated in different modes, but always adopting the 
thin 1 filter: pn in Small Window, MOS 1 in Imaging Prime Partial 
Window2 (small window), MOS2 in Timing Past Uncompressed in 
the central CCD, while in Imaging Mode in the other CCDs. The 
time resolution is 5.7 ms, 0.3 s and 1.75 ms in the pn, MOSl and 
MOS2, respectively. The data were reprocessed using version 13.5 
of the Science Analysis Software (SAS) with standard procedures. 

The RGSs were operated in spectroscopy mode and their 
data were analysed in the energy range 0.4-2.1 keV. The SAS 
task RGSPROC was used to extract the RGSl and RGS2 spectra, 
which were later combined into one single grating spectrum using 
RGSCOMBINE. 

EPIC source counts were extracted from circular regions of 
40" radius for both the pn and MOSl, adopting PATTERN from 0 
to 4 (mono- and bi- pixel events) in the pn and from 0 to 12 (up 
to 4-pixel events) in the MOSl. Background counts were obtained 
from similar sized regions offset from the source positions. Due to 
calibration issues in the MOS2 data (see Appendix lAl for details), 
we decided to use only the pn and MOS 1 data for the spectral anal¬ 
ysis, while we used all three EPIC data sets for the timing analysis. 
No further filtering based on the background level was necessary. 
Appropriate response matrices were generated using the SAS tasks 
ARFGEN and RMEGEN. 

Eree relative normalizations between the two cameras were in¬ 
cluded to account for uncertainties in instrumental responses. The 
normalization factor for the MOSl spectra relative to PN (factor 
fixed at 1) was always around 0.965±0.007. All spectral uncer¬ 
tainties and upper-limits are given at 90% confidence level for one 
interesting parameter. 

In the spe c tral fi tting we adopted the interstellar abundances 
of Wlms_et^iJj 2000l) and photoelect r ic abs orption cross-sections 


of iBalucinska-Church & McCammoij ( Il992h . using the absorption 
model PHABS in XSPEC. We checked that using a different absorp¬ 
tion model (TBabs i n XSPEC) with Wilm s et al. (2000) abundances 
and cross sections bv IVerner et al.l (Il996h similar spectral fit results 
were obtained. 


3 ANALYSIS AND RESULTS 
3.1 Timing Analysis 

The data from the three EPIC cameras in the 0.15-12 keV range 
were used for the timing analysis after conversion of their times 
of arrival to the Solar System barycenter. Using the Zl test, the 
spin period of RX J0059.2-7138 was measured to be Pspin = 

2.762383(5) s.^ _ _ 

Eollowing lKohno et al.l ( |200(]|) . we defined the pulsed fraction 

as (/max /min)/2/average, whete /max, /min, and /average ate 

the maximum, minimum, and average count rate, respectivily. The 
average pulsed fraction in the 2-10 keV energy range was 8.9 ± 
0.5 % and during the observation the value did not change more 
than ~20%. 

With respect to the last observ ation performed in 1993 
(MID 49119-49120) with ASCA, when iKohno et al.l ( l2000h found 
Pspin = 2.763221(4) s, our value implies an average spin-up rate 
of/^spin = -(1.27 ±0.01) X 10"^^ sfi-^ 

In Pig. [2 we show the EPIC pn pulse profiles in different en¬ 
ergy ranges. The pulse shape is clearly energy-dependent, with a 
broader single asymmetric peak in the 4-12 keV range, whereas 
a clear secondary pulse emerges at intermediate energies (0.5-4 
keV). Below 0.5 keV the pulsations are less evident, but still signif¬ 
icant (at 3.6 a confidence level). The pulse profile evolution with 
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Table 1. Summary of the XMM-Newton observation (ID. 0674730201) of RX J0059.2—7138 


Instrument 

Filter 

Mode 

Start and End Times 

Net Exp. time (ks) 

Net count rate (s 

pn 

Thin 1 

Small Window 2014-03-14 19:45:57 

2014-03-15 01:10:34 

13.7 

19.24±0.04 

MOSl 

Thin 1 

Small Window 2014-03-14 19:40:09 

2014-03-15 01:07:44 

19.6 

5.92±0.02 

MOS2 

Thin 1 

Timing 

2014-03-14 19:40:48 

2014-03-15 01:03:30 

19.4 

5.86±0.02 

RGSl 

- 

Specroscopy 

2014-03-14 19:39:58 

2014-03-15 01:11:55 

19.9 

0.246±0.004 

RGS2 

- 

Specroscopy 

2014-03-14 19:40:06 

2014-03-15 01:11:53 

19.8 

0.299±0.004 



0 0.5 1 1.5 


Phase 

Figure 2. RXJ0059.2-7138 pulse profiles observed with the EPIC pn cam¬ 
era in four energy bands. The vertical dashed lines and the numbers mark 
the three spin-phase intervals adopted for the phase-selected spectroscopy 
(Sect.[3H). 


energy is clearly shown also by the phase-energy image (Fig.[3l. 
This image represents the number of source counts in small phase 
and energy int ervals, normalized w ith respect to the phase-averaged 
spectrum (see lXiengo et al.|[2013h . A deficit of counts (marked by 
the white box) over a narrow range of energies (~7-8 keV) in the 
phase interval A())=0.7-0.8 is evident. This possible evidence for a 
phase-dependent absorption line will be discussed in Sect. l3.2l 


3.2 EPIC Spectroscopy 

3.2.1 Time-averaged spectrum 

We searched for possible spectral variability within the XMM- 
Newton exposure, by extracting light curves in different energy 
ranges. We found no strong evidence for intensity and spectral vari¬ 
ability on timescales longer than the spin periodicity, so we ex- 


Table 2. Results of the simultaneous EPIC pn and MOSl fit of the time- 
averaged spectram (0.2-12 keV). The double-component continuum con¬ 
sists of a cut-off power law (CUTOFFPL in XSPEC) and a black-body model. 
Both components are equally absorbed with column density Njj. Three 
Gaussian lines in emission are needed to account for positive residuals in 
the spectram (the line at E;i„ei has been added only in the pn spectrum, 
being of instrumental origin). The parameters of the black-body component 
are its temperature, kTuh, a nd it radius, ri,b. A distance of 61 kpc is as¬ 
sumed jHilditch et alj|2005l) . The cut-off power-law continuum is defined 
as A(E)=kE“'" exp — (E/Ec), where E is the photon index and Ec is the 
e-folding energy of the exponential rolloff (in units of keV). Flux;i,je^jv 
is the total flux in the Gaussian line #N. EWii„e#JV is the equivalent width 
of the line #N (eV). The unabsorbed flux is in units of erg cm“^ s“^, the 
luminosity in units of erg s“^. Li,b/L(.pi is the fraction of X-ray emission 
(corrected for the absorption) contributed by the soft component with re¬ 
spect to the cut-off power-law, in the broad energy band 0.01-10 keV. 


Parameter 

Value 

Nh (10^° cm-2) 

to 

CO 

1 + 

o o 

at b> 

kTbb (keV) 
rbb (km) 

0.093lHo5 

35011“ 

r 

Ec (keV) 

0.44l°;°2 

36l« 

(keV) 

2.23 ±0.01 

^line 2 (keV) 

^line 2 (keV) 

FluxHne2 (ph cm'^ S“^) 
EWane 2 (eV) 

0-00_o.o6 

(6±2)xl0-® 
30± 10 

EhticS (keV) 

^lineZ (keV) 

FluXHneS (ph cm~2 s“^) 
EWa„e3 (eV) 

U.yD_Q Q3 

u.uyo_o.o24 

1514 X 10-5 

30 ±8 

Lbb/Lcpl (0.01-10 keV) 

1.7% 

Unabsorbed flux (0.5-10 keV) 

1.58 Xl0-1“ 

Luminosity (0.5-10 keV) 

7 X 1037 

xj/dof 

1.405/546 


traded a time-averaged spectrum for the whole pn and MOS1 ex¬ 
posures. 

The main spectral structures in addition to a simple absorbed 
power-law component, can be seen in Fig. |4] (upper panel)’. l)-a 
soft excess below ~0.5 keV; 2)-a high energy cut-off; 3)-a narrow 
feature at 2.2 keV (present only in the pn data, very likely due to 
residual calibratio n uncertainties around the Au edge); as done by 
other authors (e.g. lPiaz Trigo et al.l20l4) we model it with a Gaus- 
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Figure 3. RX J0059.2—7138 phase-energy image (normalized to the phase-averaged spectrum) obtained from the EPIC pn and M0S2 data (the MOSl time 
resolution is too poor for such a fine phase binning). The small white box marks a hole in the image, indicating a deficit of counts, with respect to the average 
spectrum, in a narrow range of energies (~7-8 keV) and spin phases (A<?i=0.7-0.8). 


sian line, but we will not discuss it further; 4)-an emission line at 
6.5-6.8 keV; 5)-a broad positive residual around 0.9-1.0 keV. The 
presence of the latter structure also in the RGS spectrum (see be¬ 
low, Sect. l3.3l l. excludes a calibration/instrumental origin, thus we 
decided to include an additional broad Gaussian line to fit it. 

To account for these structures, we adopted a double¬ 
component continuum consisting of a cut-off power-law plus a 
black-body model together with the three Gaussian lines. The re¬ 
sults are reported in Tablej^and displayed in ¥ig.^{lower panel). 

We also tried a hot thermal plasma model (MEKAL in XSPEC) 
instead of the black-body to fit the soft excess, with similar re¬ 
sults: all the parameters were consistent with those reported in Ta- 
ble|2 except for the column density Nh=(4.1 ± 0.6) x 10^° cm“^ 
and the fraction of X-ray emission contributed by the soft com¬ 
ponent with respect to the cut-off power-law, which was 7% (in¬ 
stead of 1.7%), in the energy band 0.01-10 keV. The parameters 
of the MEKAL model were the following: a plasma temperature 
A:Tm=0.21±0.03 keV, a very low metal abundance Z < 0.007 Zq 
(too low for the SMC, see below), a normalization, normyi, of 
0.025^0 ggg cm~® (in our case, normM=2.4x 10“®^ J UeURf^V, 
where rte and nn are the electron and hydrogen densities, V is the 
emitting volume). 

We also tried to adopt two absorption models, one to account 
for the total Galactic absorption in the SMC direction (Nh fixed at 
6x10^° cm“^), the other one to account for the local absorption 
in the SMC (VPHABS model in XSPEC with an abundance fixed 
at 0.2 for metals, appropriate for the SMC), but we did not obtain 
better fits than using a single absorption model with a free column 
density. So, we decided to adopt a single model for the absorption 
(PHABS model in XSPEC). We also tried to fit the data with wa rm 
and ionized absorbers (like ABSORI in XSPECi lPone et al .Il992h in 


addition to a neutral one accounting for the interstellar absorption, 
but this did not result in better fits. 


3.2.2 Spin-phase-selected spectroscopy 

Given the energy-dependence observed in the folded spin profiles 
(Figs.|2and[3, it is not surprising that the fits to the time-averaged 
spectrum were not formally acceptable. To quantitatively explore 
the observed spectral evolution, we extracted three phase-selected 
spectra (marked by numbers in Fig. Hi for EPIC MOSl and the 
pn: A()>=0.15-0.45 (spectrum n. 1), A()>=0.45-0.85 (spectrum n. 2), 
A()i=0.85-0.15 (spectrum n. 3). 

The MOSl and pn spectra were fitted simultaneously in the 
0.2 to 12 keV range with the same model used for the time-averaged 
emission. Good fits were obtained with the spectral parameters 
listed in Table [3 Similar values were obtained using a MEKAL 
model instead of the black-body, but again with a very low metal- 
licity {Z < 0.01 Zq). 

The results reported in Table indicate that the phase- 
dependent spectral variability is mainly due to variations in the 
high energy cut-off. The iron lines around ~6.6 keV are signif¬ 
icant at 2.9 cr, 4.4 a and 4.1 a, in spectra 1, 2 and 3, respec¬ 
tively. The significance was estimated by increasing the Ax^ value 
until the confidence region boundary of the Gaussian normalisa¬ 
tion crosses 0. The iron line centroid energy shows variations from 
6.77 keV (in spec. 1) to 6.54 keV (in spec. 2). In spec. 2 the cen- 
troid energy is compatible with transitions in Fe XXI-Fe XXII ions 
jKallman & McCravlll982h . in spec. 3 with Fe XXI-Fe XXV ions, 
while in spec. 1 (line centroid at 6.77 keV), it is not consistent 
(within the 90% errors) with any known line. Although the emis¬ 
sion line in spec. 1 is significant at less than 3 ct, if we identify 
it with emission from Fe XXV at 6.7 keV, a hint for a blue-shift 
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Figure 4. RXJ0059.2-7138 EPIC time-averaged spectrum. Upper panel 
shows the fit of the pn data with a single power-law, to clearly show the 
main structures in the spectrum. Lower panel displays the photon spectrum 
from the EPIC pn and MOSl observations fitted with the model reported in 
Table[2] 


of 900 km is present, possibly indicating an outflowing pho- 
toionized plasma. Also the flux of the hroad feature at 0.9-1 keV is 
variable along the spin profile, while the other spectral parameters 
are constant, within the uncertainties. 

The phase-energy image (Fig. shows a hint for a deficit 
of counts at -^7-8 keV, with respect to the spin-averaged spec¬ 
trum. This deficit, possibly indicating an absorption feature, oc¬ 
curs in a narrow range of spin phases ((()~0.7-0.8). We extracted 
a spectrum corresponding to this phase interval and we fitted it 
adding a Gaussian line in absorption. The new fit gave a reduced 
xl=l-0l4 for 464 dof (compared to xS=l-049 for 467 dof with 
no line). The resulting line parameters were the following: cen¬ 
troid energy Eii„e4=7.50 ± 0.09 keV, width (T;i„e 4 < 0.27 keV, 
total Fluxiine 4 =—(1.281q 7 g)xl 0 “'‘ photons cm~^ s“^, equiva¬ 
lent width EWKne 4 ~—80 eV. The single trial significance of this 
line is 4.4 a, hut given that many phase intervals were explored, the 
overall significance is too small to claim a firm detection. 


3.3 RGS Spectroscopy 

The RGS first order spectra (0.4-2.1 keV) were extracted from the 
whole exposure and combined into a single grating spectrum. We 
fitted it with the same douhle-component model resulting from the 
time-averaged spectrum (black-hody plus cut-off power-law). The 
inspection of the residuals clearly showed the presence of narrow 
lines (both in absorption and in emission), together with a broad 
structure around 0.9-1.0 keV (which was also clearly present in the 
EPIC spectra). The addition of Gaussian lines to account for them, 
resulted in the line parameters reported in Table |4] The spectrum 
is shown in Fig.|5] The two brightest narrow lines, at 0.5 and 0.65 
keV, are firmly associated with N VII and O VIII Lya lines, re¬ 
spectively, while the broad component at 0.92 keV is possibly due 
to different contributions: a blend of several emission lines from 
iron in a range of ionizations states, from Fe XVIII to Fe XX, a ra¬ 
diative recombination continuum (RRC) from O VIII, likely with 
superposed weak emission lines from ionized neon (Table |4ll. The 
parameters of this broad component observed in the RGS spectrum 
are consistent with what found with EPIC. 

We then used a thermal plasma model (MEKAL) with free 
abundance instead of the black-body to account for the soft excess 
(with the metal abundance let free to vary, and resulting in a kT 
of ~0.23 keV), but in this case only the narrow line at 0.65 keV 
could be well fitted. Also adopting additional MEKAL models with 
different temperatures or a CEMEKL model (a multi-temperature 
plasma emission model built from the mekal code, with emission 
measures following a power-law in temperature) we could not ob¬ 
tain acceptable fits to the other emission lines, likely indicative of 
a photoionized nature of the emitting matter. Two absorption lines 
were also found, with possible identifications listed in Table|4] 

The O VII Hea triplet is undetected, preventing us from us¬ 
ing the so-called “G ratio” in He-like ions to obtain information 
on th e ionization process and plasma density dPorguet & DubaJ 
bOOOh . We could place the following upper limits (90% confi¬ 
dence level) to the fluxes of the three O VII triplet compo¬ 
nents: F<3.1xl0“® ph (forbidden line at 0.561 keV), 

F<5.0xl0“® ph (intercombination line at 0.569 keV), 

F<3.4xl0“® ph (resonance line at 0.574 keV). 

4 DISCUSSION 

The Be transient RX 10059.2-7138 in the SMC attracted attention 
because, like a few other X-ray pulsars, it sho wed a prominent sof t 
excess, the origin of which is still uncertain dHickox et al.ll2003l . 
Two were the main alternatives discussed by Hickox et al. to ex¬ 
plain the soft excess in RXJ0059.2-7138: either reprocessing of 
X-rays in the inner region of an accretion disc or optically thin, hot 
thermal emission. 

Thanks to the large collecting area of XMM-Newton, 
the observation obtained during the recent new outburst of 
RX J0059.2-7138, allowed us to perform a more detailed study of 
this source, even if its luminosity of Lx=7xl0®^ erg s“^ (0.5- 
10 keV) was a factor ~3 s maller than that o bserved in the 1993 
ROSAT and ASCA data dKohno et al.ll200Cih . 

The new period value, Pspin=2.762383(5) s, implies an aver¬ 
age spin-up rate of Pspin = —(1.27 ± 0.01) x 10“^^ s in the 
~20 years between the two outbursts. This spin-up value is prob¬ 
ably not representative of the true torques currently acting on the 
neutron star, since during the long time interval between the two 
outbursts, the pulsar could have been spinning-down, at least for 
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Table 3. Results of the spin-phase selected spectroscopy (spectra n. 1, 2 and 3 indicate spin phase intervals 0.15-0.45, 0.45-0.85, 0.85-0.15 respectively) fitting 
together the EPIC pn and MOSl observations in the energy band (0.2-12 keV). The model is the same adopted in Table[2] 


Parameter 

1 

2 

3 

Nh (lO^O cm-2) 


Q 1 +0.9 
^•-*^-0.8 

1<1.0 

kTbb (keV) 
rbb (km) 

0.09l°;°i 

4ZU_i3o 

0.099l°;““® 

3001““ 

0.09l°;°l 

350l=®° 

r 

Ec (keV) 

U-^o_Q Q5 

>45 

0.421°;°® 

20tt 

Q oy+0.04 
_o.04 

22l® 

FHnel (keV) 

ry qo-|- 0.03 
z.zo_0 03 

Q QQ+0.02 

Z.ZO_o 02 

9 9i::+ 0.03 
^.zo_o 03 

FHne 2 (keV) 

^line2 (kcV) 

Flux/i^e 2 (phcm“^ s~l) 
EWH^e 2 (eV) 

c ■77--I-0.06 
' -0.04 
<1.6 

(4.3±2.4)xl0-® 
23± 13 

6 154+0-05 

<0.09 

(4.9lJ-|)xl0-® 

SO-lt 

D.oy_0.08 

0 1 2+®'^^ 
^•-L^_0.08 

(8.4l|;^)xl0-® 

4Q+20 

E^ineS (keV) 

^lineZ (keV) 

Fluxii„e3 (ph cm~2 s-l) 
EWM„e3 (eV) 

0.92t°;°^ 

(40l®“)xl0-5 

80l®“ 

0.99l°;°® 

<0.09 

(9l4)xl0-® 

18l“ 

u.yo_o.o5 

0 10+0-06 
'J--LtJ_o.o4 

(16l^'‘)xl0-® 

Q-| +23 

Lbb/Lcpi (0.01-10 keV) 

0.02 

0.017 

0.015 

Unabsorbed flux (0.5-10 keV) 

1.6 xlO-io 

1.4x10-1“ 

1.7x10-1“ 

xllAof 

1.041/517 

1.067/522 

1.033/518 


Table 4. Main lines identified in the RXJ0059.2-7138 RGS spectrum (0.4-2.1 keV). Negative fluxes and negative equivalent widths (EWs) mean that the 
Gaussian lines are in absorption. 


Energyobserved 

Ion 

-^lab 

Energyiab 

a 

Flux 

EW 

(keV) 


(A) 

(keV) 

(eV) 

(10-5 ph cm-2s-l) 

(eV) 

U.OUU _0.002 

N VII 

24.779 (Lyo) 

0.50 

2.4 1?;° 

6.8±2.7 

7 ±3 

0 6516 +0-0026 
U.ODIO _0.0016 

0 VIII 

18.967 (Lya) 

0.65 

<? 

5.7±2.0 

9 ±3 

n s 11 +0.003 
U.811 _0.007 

Fe XVII (?) 

15.2610 

0.81 

5 +'^ 

^ -2 

-2.8+i:i 

-5±2 

0 Q91 +0-026 
-0.024 

Fe XVIII-Fe XX (blended?) 

- 

- 

60 +^® 

OU _4o 

I<6 

34+^0 

-3+-14 


or RRC from 0 VIII (?)- Ne IX (?) 

13.45 

0.89; 0.91 




1 0795 +o-ouuy 
L.yj/yj -0.0056 

Fe XXII (?) 

11.490 

1.079 

<9 


11 ±6 

1114 +0-013 
-0.012 

Fe XXIII (?) 

11.019 

1.125 

Q +^1 
^ -6 

-4.4j;i;“ 



part of the time. With the reasonable assumption that an accretion 
disk is present during the obs erved outburst, we can discuss the re¬ 
sults in the framework of the lGhosh & Lam9 jl979l) theory of ac¬ 
cretion torques. For a given neutron star mass, radius, and dipolar 
magnetic field, the spin period derivative, P, is expected to scale 
as (where P is the spin period, L is the X-ray luminos¬ 

ity). Assuming M = 1.4M©, R = 10km, L = 10®® erg s“^, 
B = 10^^ G the expected spin-up rate for RX J0059.2-7138 dur¬ 
ing the XMM-Newton observation is 6.8 x 10“^^ s 
(eq. 15 in iGhosh & Lamblll979l) . As expected, this value is larger 
(by a factor of ~50) than the time averaged spin-up rate between 
1993 and 2014. 

Although the pulsed fraction observed in 2013 (8.9 ± 0.5 %) 
was smaller than in 1993 (37%), the large number of counts col¬ 


lected by the three EPIC cameras allowed us to see interesting vari¬ 
ations of the pulse profile as a function of energy that were not 
evident in the previous observations. The pulse profile is double- 
peaked at intermediate energies (0.5-4 keV) and for the first time 
pulsations were seen also below 0.5 keV. In principle, the presence 
of pulsations in the soft energy range might give some information 
on the origin of the soft excess (see below), but in practice this is 
complicated by the fact that both the soft excess and the hard power 
law contribute to the counts in this energy range and could be re¬ 
sponsible for the observed modulation. 

Our observation confirms the presence of a soft excess in 
RX 10059.2-7138. Its 0.2-12 keV spectrum could be well fitted 
by a double-component continuum with a hard cut-off power-law 
(r~0.44) together with a soft excess equally well modelled by 
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Figure 5. Results of the RGS spectroscopy. The counts spectrum is shown 
when fitting the combined RGS 1 and RGS2 spectra with a double compo¬ 
nent continuum (a black-body together with a power-law) modified by the 
absorption, including six Gaussian lines to account for the brightest lines 
(four in emission and two in absorption; see Table |4). In the upper panel 
the spectrum is in units of wavelength, while in the lower panel it is in units 
of energy. Spectra have been severely rebinned in these plots, to better show 
the residuals. 


either a soft black-body (kT-^O.! keV), or a hot thermal plasma 
model (kT'~0.2 keV). While these spectral parameters are similar 
to what found in 1993, the relative contribution of the soft com¬ 
ponent to the total luminosity was much lower in 2014: the 0.1- 
10 keV luminosities of the soft and hard spectral components were 
Lbb=10®® erg s“^ and Lcpi=7.1 x lO^’’ erg s“^ in 2014, while they 
were Lsoft=8x 10®^ erg s“^ and Lhard=l-8x 10^® erg s“^ in 1993, 
implying a soft to hard luminosity ratio of 1.5% (in 2014) instead 
of 44% (in 1993). 

We discovered several features in the RX J0059.2-7138 spec¬ 
trum: narrow lines from ionized oxygen and nitrogen (in the RGS 
data), a broad feature at 0.9-1 keV, likely produced by Fe-L shell 
emission (both in the RGS and EPIC data), and an iron emission 
line at ~ 6.6 keV. The latter is clearly evident in the time-averaged 
spectrum. There is some evidence that its centroid energy varies as 
a function of the pulsar spin phase, but the significance of the line 
is marginal (< 3cr) when it attains the highest energy. 

We detected a phase dependence of the cut-off energy (from 
Ec^20 keV to Ec>45 keV), which is the main responsible for the 


spectral variability with the spin phase. A similar behavior of high- 
energy cut-off is unusual in X-ray pul sars (an example is the low- 
mass X-ray binary pulsar 4U 1626-67. ICoburn et al.l2002l) . We ob¬ 
served also a phase dependence of the flux contributed by the broad 
emission component at 0.9-1 keV, that can be explained by a vari¬ 
able emitting area illuminated/photoionized by the pulsar beam. 

An accretion disc is a natural reprocessing site for the X-ray 
emission emitted from the accreting pulsar, especially when its lu¬ 
minosity is high (~ 10 ®® erg s“®) and photons pref erentially es¬ 
cape in a fan-beam from the accretion column fe.g. IParmar et al.l 
I l989l) . fav ouring an X-ray illu mination of the equatorial regions. 
Eollowing iHickox et al.l |2004 we can assume that black-body 
emission with luminosity Lsoft=G Lx (where Lx is the total 
X-ray luminosity) results from reprocessing by optically thick ma¬ 
terial located at a distance R and subtending a solid angle Q, with 
respect to the central X-ray source. Independent on the exact ge¬ 
ometry of the reprocessing region, the distance R can be estimated 
from the relation 

„2 ^ Lx 

dvrcrT^g 

(note that this is different from the radius (r^f,) resulting from the 
normalization of the black-body component in the spectral fit (Ta¬ 
bles [2] and [3]l). If the reprocessing region is a shell at the inner 
edge of the disc, R should be of the order of the magnetospheric 
radius, Rm (or of the corotation radius if the pulsar is rotating 
close to the equilibrium period). Adopting the RX 10059.2-7138 
phase-averaged luminosity Lx=7xl0®^ erg s“^ (0.5-10 keV) and 
the black-body temperature, Tbb, of 0.09 keV (Table [tJi, we ob¬ 
tain a distance 7Z = 3x 10® cm, which is indeed very similar to the 
corotation radi us f(cor=3.3x 10® cm. A n estimate of the magneto- 
spheric radius dPavies & Pringle! 1 198 ih can be obtained from the 
relation 

Rm ~ 1.5 X 

where mi is the neutron star (NS) mass in units of solar masses. 
Re is the NS radius in units of 10® cm, L 37 is the X-ray lumi¬ 
nosity in units of 10®^ erg s“®, B 12 is the NS magnetic field in 
units 10®^ Gauss. If for RX 10059.2-7138 we assume i 3 i 2 =l we 
obtain 77^^10® cm. If the soft excess is caused by reprocess¬ 
ing of the X-ray pulsar emission by the inner edge of the accre¬ 
tion disk, some pulsations in the soft part of the spectrum are ex¬ 
pected, as indeed it was found. In this scenario, the narrow lines 
we discovered at low energies in the RGS spectrum and the ion¬ 
ized iron line at 6.6 keV can be produced by a plasma located 
in the disc atmosphere. They require the presence of matter with 
very different ionizations stages. Indeed, while plasma with highly 
ionized oxygen and nitrogen is characterized by a ionization pa¬ 
rameter ^~100 erg cm {^=LlnR?, where L is the luminosity 
of the ionizing source, n is the plasma density, R is the distance 
of the photoionized matter from the ionizing source), the pres- 
ence of highly ionized iro n Ee XXV implies ^=1000 erg cm s“® 
iKallman & McCra^ll982h . The widths we derived for the N VII 
and O VIII Lya lines are larger than the RGS energy resolution and 
imply velocities of ~ 1,000-2,500 km for the emitting plasma, 
as expected for matter rotating at distances of a few x 10 ® cm from 
the neutron star. Given the high X-ray luminosity, the presence of 
the accretion disc, and the detection of narrow lines from matter 
with very different ionization stages not compatible with a single 
plasma temperature, we can interpret the soft excess as the result of 
reprocessing from the inner edge of the disc and the narrow lines 
as due to photoionized plasma in regions above the disc. 
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An alternative interpretation is that the soft excess in 
RX J0Q59.2-7138 h entirely emitted by a thermal hot plasma 
jHickox et aljbooill . Using a MEKAL model to fit the RGS spec¬ 
trum, we could account well only for the O VIII Lya emission line, 
which we interpret as a further indication of the presence of pho¬ 
toionization, and obtained a metal abundance {Z < 0.01 Zq) too 
low for the SMC. The emission measure obtained from the nor¬ 
malization of the MEKAL model is rdV ~ 10®° cm“°. Since the 
gas density in this case should be n < 10^^ cm“°, the emitting 
spherical region for the optically thin plasma should have a radius 
f?cioud > 6 X 10^^ cm, and no pulsations in the soft excess are 
expe cted. _ _ 

iKohno et al.l ( l2000h found wavy structures below 1 keV in 
their ROSAT plus ASCA spectrum and explained them with the 
presence of a deep edge-like structure at 0.5 keV, caused by an 
over-abundance of neutral oxygen (about 6 times solar) in the cir- 
cumstellar matter. However, this was at odds with the very low 
metal abundance (0.02 solar) in their best-fit thermal plasma model. 
These authors concluded that the oxygen over-abundance is real, 
while the MEKAL model is un-physical. Also oui XMM-Newton 
spectrum showed wavy structures below 1 keV, which we fitted 
with a broad Gaussian line at ^0.96 keV, but we found no evi¬ 
dence for a neutral oxygen edge at 0.54 keV deeper than the one 
ascribable to the best-fitting Nh. We note that acceptable fits to 
the EPIC spin-phase selected spectra adopting a single absorption 
model with variable metal abundances (VPHABS in XSPEC) could 
be obtained only deleting the Gaussian component at 0.96 keV. 
However, in this case, an oxygen overabundance around 30 times 
solar was obtained, which we think is unlikely. Moreover, the re¬ 
sulting column density was very low (~5xl0^° cm~^), one order 
of magnitude lower than the Galactic one. If, instead, we use two 
absorption models, one for the Galactic absorption (PHABS) and 
one with variable abundances for the SMC absorption (VPHABS), 
the column density of the latter becomes compatible with zero and 
the fit is insensitive to the metal abundances. Given that the neu¬ 
tral oxygen edge falls within two bright emission lines in the RGS 
spectrum (N VII and O VI II), we believe that the edge-like struc¬ 
ture at 0.5 keV observed bv iKohno et al.l ( I 2 OO 0 I) was due to the low 
spectral resolution of the ROSAT plus ASCA spectra and to the 
presence of the nearby emission component at 0.9-1 keV. 

A hint of an absorption line at 7.5 keV was present only in 
a narrow range of spin phases (A((>=0.7-0.8). Since its energy is 
not consistent with any known line feature (except for the Kq, line 
from neutral Ni at 7.48 keV, which we consider very unlikely, given 
the highly ionizing flux from the X-ray pulsar), it is tempting to 
associate it with a cyclotron resonant scattering feature. In this case, 
its energy would imply a magnetic field in the scattering region of 
6.5x(l-l-z)x 10*^^ G, where z is the gravitational redshift. 


5 CONCLUSIONS 

The main results of our XMM-Newton observation obtained 
during the second bright outburst from the SMC Be/XRT 
RX J0059.2-7138 can be summarized as follows : 

• The continuum X-ray spectrum is complex and requires sev¬ 
eral components: a hard power law with a phase-dependent high- 
energy cut-off, a soft excess below 0.5 keV well fit by a black-body 
model, and a broad emission feature at ~0.9-l keV. 

• We discovered narrow emission and absorption lines from 
oxygen, nitrogen and iron in different ionization stages. 


• The relative contribution of the soft excess component to the 
total luminosity was smaller than in the 1993 observations (1.5% 
wrt 44%), when the source was a factor ~3 brighter. 

• Contrary to what observed in 1993, the pulsar pulse profile 
was strongly energy-dependent and pulsations were detected also 
below 0.5 keV. The pulsar spin period of Pspin=2.762383(5) s im¬ 
plies an average spin-up rate of Pspin = —(1.27± 0.01) x 10“*^^ s 

in the ~20 years elapsed after the only other period measure¬ 
ment of this source. 

These findings lead us to interpret the soft X-ray excess in this 
source as due to reprocessing of the pulsar emission from (part of) 
the inner edge of the accretion disc at a radius of ~ 3 x 10® cm. 
The disc inner boundary is consistent with a magnetic field of the 
order of 10^^ G which could also explain the possible absorption 
feature at 7.5 keV as an electron cyclotron line. The observed nar¬ 
row lines are most likely produced by photoionization of plasma 
located above the inner regions of the accretion disc. 
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APPENDIX A: MOS2 TIMING MODE 

Here we discuss some calibration problems that affected the 
MOS2 spectrum, which was not considered in the spectroscopy 
of RXJ0059.2-7138. Note however that these issues have no im¬ 
pact on the timing analysis, which was indeed performed using also 
MOS2 events. 

Source counts in MOS2 timing mode were initially ex¬ 
tracted from RAWX=290 to RAWX=321, resulting in a count rate 
~5 count s“^. Background was taken from an outer CCD that 
collected data in imaging mode. When compared with the spec¬ 
tra extracted from the other cameras, MOS2 spectrum showed 
a significant deficit of counts above 7 keV, a shortage becom¬ 
ing much less evident when extracting counts from a box excis¬ 
ing the central columns of the PSP (that peaks at RAWX=305). 
We extracted different spectra from boxes with a variable num¬ 
ber of excluded inner columns. Por each box region, we also 
tried different PATTERN (single-pixel events or single and dou¬ 
ble events), but this led to a marginal difference in the result¬ 
ing spectral shape in the energy range 0.6-10 keV. We fol¬ 
lowed the procedure recommended on the XMM-Newton SAS 
User Guide (issue 10.5, updated to 2014, Pebruary 7) available 
at http://XMM2.ESAC.ESA.INT/DOCS/DOCUMENTS/ to properly 
calculate the effective area associated with a spectrum extracted 
from the box in MOS2, from which the core columns have been ex¬ 
cised. Comparing these several different MOS2 spectra, and fitting 
them with a simple power law model, we found significantly differ¬ 
ent slopes, even when restricting the spectral analysis to the safer 
energy range 2-7 keV. The reason for this behaviour is unclear (we 
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did not find anything similar reported in the SAS User Guide). In 
principle, the count rate is too low (~5 count s“^) to cause possible 
pile-up problems in MOS2, although the use of the SAS task EPAT- 
PLOT on MOS2 events (for both the whole spatial RAWX region 
and the different boxes chosen) showed a significant discrepancy 
in the distributions of observed and expected event shapes in the 
energy range 0.5-2 keV, while apparently no problems were found 
above 2 keV. We also noted a significant electronic noise in MOS2 
image (TIME versus RAWX), with a reapeting pattern in RAWX 
(several equally spaced “hot RAWX columns”), appearing only be¬ 
low 0.5 keV. It is possible that this huge noise creates significant 
problems in the MOS2 spectral reconstruction. 
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